NON-LINEAR  REGRESSION  ANALYSIS:  A  GENERAL  PROGRAM 
FOR  DATA  MODELING  USING  PERSONAL  MICROCOMPUTERS 


D.  P.  NELSON 
L.  0.  HOMER 


REPORT  NO.  83-5 


s 


DTIC 

ELECTE 
APR  151983 


NAVAL  HEALTH  RESEARCH  CENTER 

P.O.  BOX  05122 
SAK  DIEGO,  CALIFORNIA  W138 


NAVAL  MEDICAL  RESEARCH  AND  DEVELOPMENT  COMMAND 

BETHESDA,  MARYLAND 


88  04  15  009 


PlBTRIBtmON  fTATEMENTT 
ApptQVMi  foi  public  icUok*: 


Afloe^lon  for _ 

¥fTS~GRA4l 
DTIC  TAB  □ 

Unannouiiood  □ 

JuBtlflontlon - 


By - 

Dlntributlon/ 

Availability  Codoa 
lAvull  Qud/or 

Blat  I  f'lpi'Liial 


NON-LINEAH  REGRESSION  ANALYSIS! 

A  GENERAL  PROGRAM  FOR  DATA  MODELING 
USING  PERSONAL  MICROCOMPUTERS 


DENNIS  P.  NELSON,  Ph.D.  * 
LCDR,  MSC,  U3N 

NAVAL  HEALTH  RESEARCH  CENTER 
P.  0.  Box  851:^2 
SAN  DIEGO,  CA  92138 


and 


LOUIS  D.  HOMER,  M.D.,  Ph.D. 

NAVAL  MEDICAL  RESEARCH  INSTITUTE 
HETHliSDA,  MD  20814 


Report  No.  03-5,  supported  by  Naval  Medical  Research  and 
Development  Command,  Betheada,  Maryland,  Department  of  the 
Navy,  under  research  Work  Unit  3M162770A871  .AD.  306  .  The 
views  presented  In  this  paper  are  those  of  the  authors.  No 
endorsement  by  tho  Department,  of  the  Navy  has  been  given  or 
should  be  inferred. 


*  Biological  Sclencos  Department 


DBnraiBunoN  statement  a 


Appiovvd  fot  public  leleoMi 
Distribution  Unlimited 


ACKNOWLEDGEMENT 

We  wish  to  thank  Dr.  William  R.  Griswold,  Department  of 
Pediatrics,  UCSD  Medical  School  for  providing  the 
experimental  data  oet  upon  which  the  examplea  are  baaed. 


SUMMARY 


\ 


\ 

This  ropot  t  docuiwenta  a  ganoral  non-linear  reg.reaalon  program  for  fitting 
data  to  non-llnaar  modeis.  Tha  program  is  basad  on  an  algorithm -by.  JiaxauaaAfe  P 
which  ua«a  a  l«a»t  aquac«»  critarlon  to  calculata  aucceaalv®  iroprovomanta  to  mn 
initial  aat  of  paramatar  eatintatoo.  The  program  ia  written  in  the  PA8IC  language 
contmon  to  mont  microcomputera^  becauae  it  la  ea»y  to  use  and  to  transport  between 
machines  from  different  manufacturers.  The  majority  of  InexponBive 
microcomputera  do  not  offer  matrix  operations  as  part  of  their  BASIC  Interpreter. 


The  program  preaented  hare,  therefore,  suppllee  aubroutinoa  in  BASIC  for  tha 
lerolng,  tranapoaing  and  itivertlng  of  the  required  raatricea,  to  make  it 
compatible  with  moat  microcomputera  available  today.  The  report  givoa  oxamplea 
and  program  output  beeed  on  a  d««monal.ration  data  set  In^'olving  antigan—antibcAy 
oomploxation  in  uolutiim.  Two  derivr.tions  of  function  eubroutinaa  are  given  to 
aaslMt  the  uaer  in  developing  his  own  functiona.  A  complete  liatlrig  of  the 
neoesaary  progtama  le  glvan  along  with  a  aectlon  on  program  cautlona. 


INTRODUCTION 


The  advent  of  the  pernonal  computer  and  the  development  of  general  purpooe 
software  has  now  made  It  poaalble  for  blo-medlcal  scientists  to  take  advantage  of 
powerful  modeling  tools  previously  available  only  to  mathematicians  or  computar 
ocientists.  Care  must  be  taken,  however,  to  insure  that  programs  are  correctly 
applied  and  that  any  limitations  in  the  methods  are  clearly  understood. 

One  modeling  application  which  appears  repeatedly  in  bio-medloal  research  is 
that  of  predicting  an  outcome  on  the  baais  of  experience.  This  statlstloal 
method,  known  as  regression  analysis,  requires  that  a  functional  relationship 
between  the  dependent  and  independent  variables  be  specified.  In  the  past, 
regreeeion  analysis  has  been  largely  limited  to  linear  models.  Such  models  can 
be  aolved  by  hand  with  a  single  matrix  inveralon,  although  they  are  more  eaaily 
solved  using  a  computer.  Certain  other  models  can  be  made  linear  by  parameter 
transformation  in  order  to  utilise  linear  regraeslon  techniquea.  Tranaformatlons 
can  introduce  unwanted  and  sometimes  unsuspected  limitations  or  assumptions  Into 
the  model  and  must  be  used  with  care.  The  majority  of  models  encountered  In 
blo-msdical  research  and  In  problems  of  solution  equilibrium,  however,  are 
non-linear  and  many  cannot  be  transformed  into  linear  form.  In  order  to 
adequately  model  these  systems  a  requirement  exists  for  a  general  purpose 
non-linear  regression  technique  which  la  both  eufficlently  general  and  robust  to 
be  applicable  to  a  wide  variety  of  research  intereata. 

This  report  documents  a  general  non-linear  regression  program  for  fitting 
data  to  non-linear  models.  The  program  ie  written  In  the  BASIC  language  common  to 
most  microcomputers.  The  program  has  been  adapted  for  use  on  several, 
stand-alone,  microcomputer  systems  including  the  Wang  2200  aerlas,  tha  Radio 
Shack  TRS-60  Modela  I  and  II,  and  the  Apple  II.  Because  SASIC  la  common  to  most 
microcomputer  systems  the  program  Is  transportable  between  machines  with  relative 
ease.  Some  microcomputers  offer  matrix  operations  as  part  of  their  BASIC 
interpreter,  however  the  majority  of  inexpensive  microcomputers  do  not.  The 
program  preaentad  here  supplies  subroutines  in  BASIC  for  the  raroing,  transposing 
and  inverting  of  the  required  matrices,  to  make  it  compatible  with  moat 
mlorocomputara  available  today.  If  the  program  ia  to  ba  used  on  equipment  which 
supports  matrix  operations,  the  matrix  subroutine  calls  can  he  replaced  with  the 
correepondlng  matrix  keywords. 


METHOD 


The  program  ie  baaad  on  an  algorithm  by  Marguardt  (1)  which  uses  a  Taylor's 
series  expansion  to  give  successive  improvements  to  an  initial  set  of  parameter 
estimates.  The  method  is  actually  a  compromise  between  the  Taylor's  series 
(linearization  method!  and  the  method  of  steepest  descent  (gradient  method).  It 
combines  the  best  features  of  both  methods  while  avoiding  their  most  serious 
limitations  (2).  It  shares  with  the  gradient  methods  thoir  ability  to  converge 
from  an  initial  guess  which  may  be  outside  the  region  of  convergence  of  other 
methods,  and  with  the  Taylor  series  methods  their  ability  to  rapidly  converge 
once  in  the  vicinity  of  the  minimum,  An  attenuation  parameter  (lambda)  is  used 
to  interpolate  between  the  two  methods  as  needed.  Making  lambda  large  favors  the 
gradient  method,  and  expands  the  region  of  convergence;  malting  lambda  small 
selects  the  Taylor  series  and  favors  rapid  convergence.  Although  no  single 
method  can  be  considered  best  for  all  non-linear  problems  Marquart's  method  is  a 
sensibLe  first  choice,  provided  the  Initial  parameter  estimates  are  reasonable, 

The  user  supplies  a  subroutine  which  specifies  the  calculation  of  an  error 
between  an  observed  variable  and  a  model  of  the  observations,  The  model  is  a 
mathematical  relationship  between  the  independent  variables  (X(J,1)-  XlJ,5))>  the 
coefficients  or  parameters  being  estimated  ( B  (1 ) -B(  10 )  ) ;  and  the  independent 
variable  (Y(J)).  The  subscript  (J)  is  used  by  the  program  to  indicate  each 
distinct  data  point  and  must  be  Included  in  the  model.  The  user  provides  the 
data  and  the  starting  guesses  for  the  parameters  or  coefficients  (B's)  being 
estimated  by  t'.ie  least  squares  criterion.  The  program  then  calculates  a  sum  of 
squared  errors  (SSE).  The  program  also  calculates  numerical  partial  derivatives 
of  the  function  with  respect  to  the  B's  and  uses  these  together  with  the 
attenuation  parameter  (lambda)  to  select  a  new  set  of  B's.  If  the  new  B's  give  a 
better  SSE  then  lambda  is  divided  by  10  and  the  cycle  is  repeated.  If  the  new 
B's  do  not  Improve  the  SSE,  then  lambda  is  multiplied  by  10  and  another  set  of 
B's  is  calculated.  At  each  Iteration  the  new  parameters,  the  SSE,  lambda,  and 
R-squared  are  printed  out.  All  this  is  repeated  until  the  pre-selected  number  of 
cycles  has  been  exhausted  or  until  auccassive  changes  in  the  B's  are  less  than 
one  part  in  10,000  (Line  470).  Detailed  accounta  of  this  method  may  be  obtained 
from  the  references  listed. 


THE  PROGRAM 


The  progcam  actually  exists  In  fout  Baparats  p'v.'rtu  (listed  In  Appendix  A) 
which  are  stored  as  Independent  disk  files.  As  isted  thov  run  on  a  Radio  Shack 
TRS-80  Modal  II  but  can  easily  be  adapted  to  any  machine  using  Microsoft  BASIC. 
The  program  statements  that  differ  from  machine  to  machine  ure  those  concerned 
with  I/O.  On  some  machines  the  multiple  statement  lines  will  need  to  be  listed 
as  single  statements.  Although  several  versions  of  this  program  have  contained 
subroutines  for  plotting  the  data  points  as  wall  as  the  functional  relationship 
in  order  to  visualize  the  model  “fit")  plotting  routines  have  been  omitted  from 
this  report  because  they  are  extremely  machine  dependent.  If  plotting  is  desired 
It  is  usually  quite  simple  to  add  plot  routines  for  a  specific  machine.  The 
experimental  data  points  are  contained  In  the  variable  arrays  YO  and  X(,).  A 
aeries  of  X  values  can  be  generated  incrementally  between  a  lower  and  an  upper 
limit,  and  a  corresponding  Y  value  can  be  calculated  for  each  using  the  function 
subroutine  in  Line  50  of  the  program.  The  model  curve  can  then  be  plotted  and 
compared  to  the  plot  of  the  experimental  data.  Of  course  only  one  independent 
variable  can  be  plotted  at  a  time. 

The  first  program  is  called  "NOMLIN/DAT"  and  is  used  to  facilitate  data 
entry  and  to  create  disk  data  files.  Many  Independent  data  sets  can  be  created 
and  stored  on  disk  under  different  names  using  "NONLlN/DAT" .  The  main  program  ia 
called  "NONLIN/REG"  and  performs  the  actual  regression  analysis.  The  third  la 
"NONLiM/pRT"  and  cannot  be  used  independently  buc  may  bo  merged  with  "NONLIN/REG" 
to  divert  the  output  from  the  screen  to  the  printer.  The  fourth  Is  the  function 
subroutine.  It  is  called  by  any  8  character  description  "********'•  plus  the  file 
extension  /RUN.  The  function  subroutine  Is  a  BASIC  merge  file  and  must  be  merged 
with  "NONMN/REQ"  prior  to  program  execution.  This  design  allows  many  functions 
to  be  stored  on  the  disk  with  only  a  single  copy  of  the  "NONLIN/RBQ"  program. 
Since  both  "NONLIN/PRT"  and  "********/kun"  are  merge  files  they  must  be  saved  to 
disk  in  ASCII  format. 

Program  1  "NONLIN/DAT"  and  program  3  "NONLIN/PRT"  are  self  explanatory, 
require  no  changes  and  should  be  used  as  listed  except  for  any  machine  apacifio 
modifications.  Program  2  "NONLIN/REG"  has  the  option  of  adding  a  "Plot  Data 
Points"  routine  between  lines  1006  and  1  499  and  a  "Plot  Leaat  Squares  Fit" 
between  lines  1506  and  1994.  A  subroutlna  to  scale  the  data  or  to  establish 

5 


1 


upper  and  lower  limits  for  the  plot  can  be  added  In  this  same  region.  It  ia 
possible  to  change  the  convergence  criterion  by  changing  the  value  in  line  470  to 
a  larger  or  amaller  number.  The  rest  of  the  program  is  quite  general  with  the 
ability  to  handle  up  to  100  data  points  and  up  to  5  independent  variables 
(X ( J ,  1 )-X( J, 5 ) )  automatically.  The  total  number  of  data  points  and  the  number  of 
X  variables  are  read  in  with  the  data  set.  The  number  of  parameters  la  governed 
by  the  function  and  is  limited  to  10  (B(l )-B<10 ) ) .  If  computer  RAM  storage  la 
limited  it  may  be  necessary  to  reduce  the  size  of  the  X( , )  and  Y(l  arrays  to 
conserve  space.  The  progrem  Itself  occupies  about  5400  bytes  of  RAM  and 


to  6000  bytes  when  the  printer  routine  and  an  average  function  routine 


are  added.  Program  space  can  be  reduced  to  about  4400  bytes  by  compression  if 
absolutely  necessary.  This  Is  not  recommended >  however,  beoauie  the  program  code 
becomes  unreadable  and  the  merging  of  printer  and  function  routines  impossible. 
The  data  arrays  require  an  additional  5500  bytes  for  a  total  program  raquirsmont 
of  11,500  by  tea.  The  beat  way  to  recover  space  is  to  redefine  the  data  arrays, 
it  possible,  to  use  less  memory.  The  majority  of  microcomputars  can,  however, 
accommodate  this  program  eaaliy  in  its  entirety.  In  general  programs  1-3  require 
little  operator  intervention  or  attention. 

Program  4  ******** Is  the  only  part  of  the  program  that  is  really 
variable  and  must  be  written  by  the  operator.  This  is  the  mathematical  function 
or  model  to  which  the  data  are  to  be  fitted.  The  subroutine  begins  in  line  50 
and  uan  continue  through  line  99.  It  is  located  in  the  beginning  of  the  program 
for  easy  viulbllity  and  because  it  is  the  most  frequently  called  routine  in  the 
program.  Since  the  BASIC  interpreter  starts  searching  for  subroutines  at  line  1, 
the  location  of  the  function  subroutine  early  in  the  program  provides  a  slight 
speed  advantage.  The  subroutine  is  created  like  any  other  BASIC  program  and  then 
saved  to  disk  in  ASCII  format  so  that  it  can  bs  merged  Into  the  main  program. 
Since  BASIC  doee  not  differentiate  between  variables  in  the  main  program  and 
those  in  subroutines,  care  must  be  taken  not  to  use  any  variables  in  the 
subroutine  which  ace  uned  in  the  main  program.  A  variable  cross  refersnce  index 
is  given  in  Appendix  D  to  help  users  avoid  main  program  variables  In  the  function 
subroutine.  Array,  string  and  explicitly  defined  variables  of  the  same  name, 
however,  are  all  distinct  In  Microsoft  BASIC.  Appendix  C  gives  the  derivation  of 
the  function  subroutines  which  oorrespond  to  those  used  for  the  example  program 
output  in  Appendix  B.  The  function  routine  has  only  two  tasks,  to  use  the 
parameters  (B's)  and  the  X  values  to  caloulats  a  y(aalo)  value  and  to  compare 


thlf  valu0  with  the  true  y  value  to  determine  the  error,  Thla  error  rauet  be 
Btored  In  the  variable  E  (E>y( J)-Y{calo) >  and  will  be  uaed  by  the  main  proqren  to 
compute  the  8SE,  The  subroutine  must  then  return  to  the  main  program.  It  Is  a 
good  praotlae  to  print  the  function  subroutine  (LLIST  50-$9)  before  running  any 
program.  This  provides  a  permanent  record  of  the  model  used  for  that  fit  on  the 
sane  page  as  the  program  output. 


RONNIMG  THE  PROGRAM 


Host  mlorocomputera  first  load  a  Dish  Operating  System  (DOB).  To  run  a 
BASIC  program  the  BASIC  Interpreter  must  be  loaded  before  the  program.  Allowance 
for  at  least  one  I/O  buffer  must  be  made  when  BASIC  la  loaded.  The  sequence  for 
loading  NOMLIN  ist  LOAD  “NOMLIN/RSO" )  MERGE  "NONLIN/PRT"  If  the  output  lu  to  go 
to  the  printer)  and  MERGE  ■e******VEUM''  (the  function  subroutine  whioh  has  been 
previously  created  and  stored  on  dlslt).  Examples  of  program  output  for  one,  two 
and  three  parameter  models  using  the  same  data  are  given  In  Appendix  B.  Typing 
RUN  begins  execution  of  the  program.  The  program  then  requests  the  following 
Inputs* 

Data  File  Name?  Response.i  Enter  the  name  given  the  data  file  at  the 

time  of  file  creation. 

Max.  No.  lLuidi.lona7  —  ncsponssi  Enter  a  limiting  number  of  iterations 
in  case  the  model  does  not  converger  a  good  practical  number  is  15. 

Initial  Value  Lambda?  — -  Responset  Enter  the  starting  value  for  lambdar 
usually  1.  The  more  accurately  the  atarting  B's  are  Icnown  tha 
smaller  lambda  can  be. 

No.  of  Parameters?  -  Responsai  Enter  number  of  parameters  in  the 

regression  model.  Must  be  thn  soma  as  the  f  of  parameters  used  in 
the  function  B(l)  to  D(IO). 

Starting  Parameters? - Responast  Enter  the  beat  eetimate  possible  for 

each  of  the  parametera  used  in  the  function.  Ifr  for  example,  3  was 
entered  under  No.  of  Paremeters  than  s  separate  estimate  must  be 
given  for  each  of  the  3  paremotera  here.  It  is  possible  to  defeat 
the  program  by  entering  unrealistic  values  for  the  It's  at  this  point. 
The  eatlmates  must  correspond  exactly  aa  used  In  the  modal  B(l), 
B(2)r  B(3)  eto. 

Print  Data? - Rsaponaei  Typing  1  will  oausa  the  data  set  to  be  printed 

out,  Y  first  then  each  of  the  five  X's.  If  any  of  tha  X'a  are  not 
uaad  a  0  la  prlntad  instead. 


Prom  this  point  on  tho  program  runs  automatically  printing  the  number  for 
each  Iteration,  the  values  of  the  parameters,  the  S8B,  the  R-aquare  and  lambda. 
Only  those  steps  In  which  the  8SE  improves  count  aa  an  Iteration.  If  the  BSE  has 
Inoroasad  instead  of  decreased,  a  new  set  of  B's  will  b*  found  using  the  old  B's 
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and  a  10  fold  larger  value  for  lambda.  Lambda  la  not  printed  and  the  Iteration 
counter  la  not  incremented.  Iteration  contlnuea  until  either  the  convergence 
criterion  la  met  or  the  program  runa  out  of  Iteratlona  aa  apeolfled  on  program 
entry.  It  then  seta  lambda^O  and  caloulateo  a  final  set  of  parametera/  a  final 
SSB  and  a  final  R-aquare.  It  also  gives  the  overall  variance  and  the  atanderd 
deviation  and  an  approximate  standard  error  for  each  parameter  with  its 
coefficient  of  variation.  The  more  nearly  linear  the  model,  the  more  accurate 
are  these  approxlmationa. 

The  program  then  offers  the  option  to  print  out  the  var lance-covarlanco 
matrix,  the  correlation  matrix  and  a  table  of  residuals.  The  table  of  residuals 
la  useful  for  quickly  detecting  bias  In  the  model,  and  the  correlation  matrix  for 
spotting  Interactiona  between  parameters.  The  program  alao  offeta  a  complete 
printout  of  the  X  and  Y  data  and  uaes  the  beet  values  of  the  paramoters  (B's), 
found  by  the  Iteration,  to  compute  the  corresponding  Ytcalc)  values  and  their 
approximate  standard  errors.  The  95%  confidence  Interval  on  V(oalc)  is  also 
printed  as  (4-/-  2  SD)  as  Is  the  residual.  Only  the  first  value  of  X  Is  printed. 
If  plotting  subroutines  have  been  added  to  the  program,  plots  of  the  data  and 
function  can  be  drawn  at  this  time.  The  program  then  either  proceeses  another 
data  set  or  terminates. 

EXAMPLES 

Three  separate  program  outputs  for  a  one,  a  two  and  a  three  parameter  model 
are  provided  in  Appendix  B.  The  same  data  set  Is  used  for  all  three  examples  to 
demonstrate  how  changing  the  mathematical  model  to  more  accurately  represent  the 
actual  epeclea  In  solution  results  in  a  considerable  improvement  In  the  fit  of 
the  experimental  data.  Derivation  of  the  functions  uaed  for  these  models  Is 
provided  In  Appendix  c  to  Illustrate  two  common  methods  for  developing  the  model 
functions.  The  data  set  used  for  the  examples  represents  a  number  of  serial 
dilutions  of  a  solution  containing  an  antigen-antibody  oomplex.  The  Y  value 
represents  that  fraction  of  the  total  antigen  bound  to  antibody.  The  YCJI'e 
were  obtained  by  dividing  the  amount  of  radioactivity  present  in  an  ammonium 
sulfate  precipitation  of  the  complax  by  the  total  radioactivity.  TheX(J,l) 
value  is  the  di  lution  factor  for  each  of  the  11  solutions,  The  total  starting 
concentrations  for  both  antigen  and  antibody  are  known.  The  total  concentrations 
for  each  dilution  are  obtained  by  multiplying  the  total  concantratlona  by  the 


dilution  factor  X(J,1).  You  will  note  that  the  undiluted  solution  is  33.4%  bound 
while  the  1/50  dilution  is  only  9.3%  bound. 

The  first  example  in  Appendix  B  is  the  111  model  for  anti<jen-Antlbody 
binding.  Only  one  parameter  Is  needed  for  this  model  namely  the  equilibrium 
binding  constant  K.  This  model  is  the  slmpleet  case  for  a  solution  interaction 
since  it  permits  only  a  a'ngle  species  of  both  antigen  and  antibody  and  only  one 
interaction  between  them.  Models  of  this  type  are  relatively  easily  solved 
algebraically  in  terms  of  known  quantities  (ouch  as  the  total  antigen  and  total 
antibody  concentrations)  and  the  parameters  being  estimated  (see  Appendix  C). 
The  equilibrium  expression  and  the  mass  balance  equations  reduce  to  a  quadratic 
equation  In  either  free  antigen  or  free  antibody  concentration  which  can  be 
solved  using  the  quadratic  equation  aa  is  done  hero.  This  method  is  preferable 
to  iterative  methods  for  function  eolution  since  the  evaluation  ia  direct,  not  by 
successive  approximation.  There  are  many  instancasr  particularly  with  multiple 
equilibria,  where  higher  order  equations  are  encountered  which  are  more  difficult 
to  solve  algebraically.  In  these  cases  numerical  procedures  for  evaluating  the 
function  are  preferred  although  computer  time  may  be  lengthened  substantially. 
The  numerical  solution  of  a  function  is  shown  in  example  3.  Many  problems  in 
solution  equilibrium  can  be  expressed  in  terms  of  the  fraction  of  oome  component 
bound  aa  a  function  of  an  unbound  species.  This  type  equation  often  reduces  to 
the  form  ( Y-B  *X/ (14-B*X )  where  Y  is  the  fraction  of  a  speoles  bound,  B  is  a 
parameter,  usually  a  binding  constant,  and  X  ia  the  free  concentration  of  the 
binding  species.  The  single  parameter  model  yields  a  binding  constant  of 
S.lE-tS,  with  a  standard  error  of  about  171.  The  SSB  is  .0323  and  the  R~square  of 
.43  shows  that  only  434  of  the  total  variation  is  explained  by  the  model.  The 
table  of  reelduale  ehows  a  clear  bias  in  the  fit  with  large  negative  errors  at 
one  end  of  the  data  to  large  positive  errors  at  the  other.  This  is  a  clsur 
indication  that  the  model  le  not  quite  correct. 

The  second  example  ie  an  attempt  to  improve  tho  fit  by  introducing  a  second 
parameter  but  using  the  same  lil  binding  model.  Here  we  postulate  that  the 
starting  total  antibody  concentration  in  incorrect  and  that  some  other 
"effective"  antibody  concentration  more  accurately  represents  conditions  in  the 
solution.  The  first  parameter  is  the  binding  constant  and  the  second  is  the 
"effective"  total  antibody  concentration.  As  you  can  see  addition  of  this 
parameter  reduces  the  SSE  to  .0016  and  improves  the  R-square  to  97%  explained 
variation,  a  clear  Improvement  in  the  fit.  The  binding  constant  is  changed  to 


2.7E-f6  with  a  atandard  error  of  12%.  The  total  antibody  concentration  drops  from 
4.5E-6  to  an  effective  concentration  of  2.3E-6  which  greatly  improves  the  fit  to 
the  Itl  modal.  The  standard  error  in  this  concentration  parameter  is  only  4%. 
The  residuals r  although  smaller  than  those  In  example  1,  still  indicate  a  bias 
since  they  are  first  positive,  then  negative  and  then  positive  again.  In  general 
one  hopes  for  randomly  distributed  signs  on  the  residuals.  The  correlation 
matrix  Indicates  a  fairly  high  correlation  between  the  two  variables  at  .84.  On 
statistical  grounds,  an  increase  in  the  number  of  degrees  of  freedom  (parameters) 
in  the  model  la  expected  to  decrease  the  SSE  and  Improve  R^square.  whether  or 
not  the  decrease  In  the  SSE,  on  addition  of  a  parameter,  is  significant  can  be 
determined  using  an  F-teat. 

Example  3  expands  the  model  from  a  lil  single  species  binding  model  to  a 
multiple  species  model,  This  3  parameter  model  postulates  that  the  total 
antibody  concentration  is  divided  between  two  different  antibody  species,  each 
with  a  different  binding  constant,  competing  for  a  single  species  of  free 
antigen.  The  Inverse  model  where  two  different  spscles  of  antigen  compete  for 
free  antibody  also  fits  the  data  well.  The  first  two  parameters  are  the  two 
binding  constants  and  the  third  Is  a  distribution  coefficient  which  determines 
what  fraction  of  the  total  antibody  belongs  to  each  species,  These  equilibrium 
expressions  and  the  mass  balance  equations  would  reduce  algebraically  to  a  cubic 
equation  in  the  known  guantltles  and  the  parameters.  In  general  the  resultant 
equation  la  a  polynomial  of  order  (#  binding  constants  +  1).  Solving  higher 
order  equations  in  terms  of  Icnowns  and  finding  their  roots,  Is  more  Involved  than 
the  solution  of  a  quadratic  in  examples  1  and  2.  since  equilibrium  equations  are 
positive  definite  functions,  an  easier  method  is  available,  although  it  requires 
more  computer  time,  it  is  possible  to  use  a  binary  search  to  find  the  true  value 
for  free  antibody  concentration  by  successive  approximation. 

The  binary  search  wor)(a  as  follows.  First  one  uses  the  equilibrium  and  mass 
balance  equations  to  solve  for  the  free  concentrations  of  antibody  1  (Bl)  and 
antibody  2  (B2)  in  terms  of  the  binding  constants  (Kl)  and  (K2),  the  total 
antibody  concentration  (TB)  and  the  free  antigen  concentration  (G>.  We  alao 
solve  for  total  antigen  concentration  (TG)  in  terms  of  the  (Kl,  K2,  Bl,  B2)  and 
(G).  All  except  (G)  are  known.  We  can  find  (G)  by  guessing  a  concentration  half 
way  between  zero  and  the  total  antigen  concentration  (TG).  We  then  use  this 
guess  to  calculate  a  value  for  tg,  TG(calc).  We  compare  this  with  the  true  value 
for  TG.  If  TG(calc)  1.8  too  high  then  the  guess  for  G  was  too  high  and  we  set  the 
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now  upper  limit  on  G  to  the  old  value  for  Q  and  keep  the  lower  limit,  If  TG  is 
too  lew  then  G  was  too  low  and  we  set  the  lower  limit  on  G  to  the  old  Q  value  and 
keep  the  upper  limit.  A  new  G  is  then  ohoeen  half  way  between  the  new  upper  and 
lower  limit  i  and  th9  procedure  la  repeated.  The  search  continues  until  the  new 
value  for  G  does  not  differ  signif ioantly  from  the  old  value.  The  function  then 
uses  (TG-0)/TG  to  calculate  the  fraction  antigen  bound,  Vlcalo),  computes  the 
error  E  and  returns. 

The  results  yield  a  high  affinity  antibody  with  a  constant  of  4.4E+6,  and  a 
low  affinity  species  with  a  binding  constant  of  6.0E't4.  These  constants  have 
coefficients  of  variation  of  12%  and  20%  respectively.  The  third  parameter  tells 
us  that  the  proportion  of  high  affinity  antibody  la  39%  and  low  affinity  antibody 
61%.  You  will  also  note  that  the  Sbu  has  been  decreased  to  .00033  and  the 
R-squara  increaaad  to  99.4%  explained  varl(>tion.  Thle  la  a  significant  | 

improvamant  in  fit  ovar  the  two  parameter  model  in  example  The  residuals  alao 
look  vary  good.  They  are  all  small,  have  randomly  dlsbrlbutad  signs  and  show  no 
particular  bias.  The  correlation  matrix  shows  a  correlation  of  .45  between  the 
two  binding  constants,  .04  between  Xl  and  the  distribution  constant  (R)  and  .62 
between  K2  and  R,  On  the  basis  of  the  random  residuals,  the  low  SSE  and  the  high 
R- square  we  would  conclude  that  the  model  in  example  3  adequately  represents  the 
data,  Although  there  may  be  other  models  which  also  fit  this  data,  2  antibodies 
1  antigen  is  a  good  model. 


PROGRAM  CAUTIONS 


The  Taylor's  series  used  in  the  program  has  a  region  of  convergence.  Making 
lambda  large  Increases  the  else  of  this  region,  however  it  is  still  possible  to 
give  Initial  parameter  estimates  which  ara  outside  this  region.  When  this 
happens  the  program  may  diverge  instead  of  converging  and  the  program  will 
terminate  with  an  error  condition,  either  a  division  by  sero  or  an 
overflow/underflow  error.  This  can  alao  happen  if  the  model  la  incorrectly 
postulated  or  an  algebraic  error  is  made  in  the  function  derivation.  It  can  also 
happen  if  there  is  an  error  In  the  data.  If  a  program  error  happens,  the  data 
should  be  carefully  checked  for  accuracy,  the  functional  model  checked  for 
algebra  and  the  starting  parameters  checked.  If  this  doesn't  help  a  new 
parameter  set  should  be  tried  or  a  larger  starting  lambda  value. 
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Many  functions  of  Interest  to  the  researcher  are  limited  In  their  range  by 
nature.  For  example  concentrations  and  therefore  binding  constants  can  only  have 
positive  values  (negative  concentrations  are  undefined).  The  regression  program 
has  no  such  limitation)  however)  since  mathematloally  negative  and  positive 
numbers  are  equally  valid.  This  can  present  a  problem  since  the  program  can  find 
solutions  with  negative  (meaningless)  as  well  as  positive  binding  constants.  The 
problem  can  be  avoided  by  careful  restructuring  of  the  function  1 .  e.  by 
transforming  one  or  more  parameters  to  new  parameters  with  natural  limits 
corresponding  to  the  desired  limits.  For  example  one  can  use  the  log  K  value 
Instead  of  the  K  parameter  since  the  log  function  is  naturally  llmltad  to 
positive  numbers.  In  example  3  we  used  another  technique  to  limit  the  range  of  a 
parameter.  Instead  of  using  the  parameter  R  for  thg  distribution  coefficient 
between  the  antibody  species/  ws  used  the  squared  Sin  function.  Since  the  only 
meaningful  values  for  R  are  positive  numbers  between  0  and  1  the  range  of  R  must 
be  restricted.  This  is  done  by  letting  (R'-Bin(B(3)  )*Sin(B(3) )  where  the  new 
parameter  Is  the  argument  of  the  Sin  function  In  radians.  Now  R  la  llmltad  to 
the  positive  range  ()  to  1  whereas  if  It  were  a  parameter  Itself  it  could  Include 
all  positive  and  negative  real  numbers.  The  immediate  mode  of  BASIC  can  be  used 
after  completion  of  the  program  to  print  the  actual  value  of  the  parameter  of 
Interest. 

One  must  also  recognize  that  just  because  the  program  finds  a  parameter  fit 
to  a  data  set  does  not  make  that  model  true.  The  mathematical  representation  may 
be  useful  for  prediction  but  may  or  may  not  shad  any  light  on  the  actual 
mechanism  of  the  process  under  Investigation.  The  model  should  not  only  be 
algebraically  correct  it  should  also  have  either  physical  or  chemical 
significance/  l.e.  It  should  represent  and  be  consistent  with  all  other  known 
oharacteriatics  of  the  real  system. 

Certain  models  and  data  sets  have  an  SSB  surface  which  exhibits  regional  as 
well  as  global  mlmlma.  If  this  occurs  it  Is  possible  for  the  program  to  locate  a 
local  minimum  Instead  of  the  global  minimum.  in  this  case  the  convergence 
criterion  Is  met  but  the  program  has  fallen  into  a  pit  and  has  not  located  the 
lowest  point  on  the  surface.  To  test  for  this  condition  one  should/  when 
initially  Investigating  a  new  model/  examine  the  surface  using  several  starting 
values  for  the  parameters.  If  the  program  converges  to  the  same  point  from 
several  different  atartlng  points  one  can  be  reasonably  confident  that  the  global 
minimum  has  been  found. 
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APPENDIX  A 


PROGRAM  (1)  NONLIN/DAT 

10  REM  ****  DATA  INPUT  ROUTINE  FOR  NONLIN  **** 

IS  DIM  y(100),X(100,5) 

20  CDS  I  INPUT  "NUMBER  OF  DATA  POINTS  (100  MAX)*;N 
25  INPUT  "NUMBER  OF  INDEPENDENT  (X)  VARIABLES  (5  MAX)"}Nl 
30  INPUT  "NAME  OF  DATA  FILF*;Z$ 

3b  PRINT  t  PRINT  "Y  VAR.  ,  X  VARS." 

40  PRINT  " - "  I  PRINT 

45  FOR  I-l  TO  N 

50  PRINT  "ENTER  Y<"|Ij")">  t  INPUT  Yd) 

55  FOR  J-l  TO  Nl 

60  PRINT  "X("jl|jri")"»  t  INPUT  X(I,J) 

65  NEXT  J,1 
70  OPEN  "0"il(Z$ 

75  PRINT  »1,  N,N1 
SO  FOR  I-l  TO  N 
85  PRINT  #1,  Yd) 

90  FOR  J-1  TO  Ml 
95  PRINT  tlr  Xd.J) 

100  NEXT  J.I 
105  CLOSE  1 
110  END 
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PROGRAM  ( 2 )  NONLIN/REQ 

10  CLS  I  PRINT  "****  NON-LINEAR  REGRESSION  ANALYSIS  *«**"  i  PRINT 
15  DBPINT  I,.TfK 

20  DIM  B<10),a(.lO),P(10)  ,Y(100)  ,X(100,S) 

25  DIM  A(10, 10), AA(10, 10), 0(10, 10) ,0(10,10) 

30  GOSUB  2400 
35  GOTO  110 

50  REM  ***  INSERT  FUNCTION  HERE  (LINES  50-99)  *** 

the  function  must  CALCULATE  A  VALUE  (F)  FOR  y(J), 

***  COMPUTE  THE  ERROR  BETWEEN  THE  OBSERVED  AND  CALCULATED 
***  Y  VALUES  (E-Y(J)-Pl,  THEN  RETURN. 

100  INPUT  "NEW  POINT  SET  (O-NO,  l'-YES)";Ee  :  PRINT 

105  IP  E6>0  THEN  RUN  ELSE  1995 

110  INPUT  "MAX.  NO.  OP  ITERATIONS" >T1 

120  INPUT  "INITIAL  VALUE  OP  LAMBDA"}L 

130  INPUT  "NO.  OF  PARAMETERS" ;Nl 

140  PRINT  I  PRINT  "STARTING  PARAMETERS"  i  PRINT 

145  FOR  I-l  TO  N1 

150  PRINT  «B(">I|")";  l  INPUT  B(I)  i  NEXT  I 

160  PRINT  I  INPUT  "DO  YOU  WANT  TO  PRINT  DATA  (O-NO,  1-YBS)">B5 
170  IP  E5-0  THEN  200 

175  PRINT  I  PRINT  "OBSERVATIONS  Y,X(1-S)"  l  PRINT 
180  FOR  I«1  TO  N 

190  PRINT  "("Hi")"jTAB(7)»Y(I)|TAB(17)»X(l,l )»TAB(27);X<I,2)» 
TAB(37)|X(I,3))TAB(47)rX(I,4);TAB(57)|X(I,5) 

195  NEXT  I 
200  Zl-0  !  22-0 
205  FOR  I-l  TO  N 
210  Zl-Zl+Yd) 

215  Z2-Z2'fY(I)"Y(I) 

220  NEXT  I 
225  T2-0 
230  T2-Ta+1 
235  SO-0 
240  GOSUB  2000 
245  GOSUB  2025 
250  GOSUB  2050 
255  GOSUB  2075 
260  FOR  J-1  TO  N 
265  GOSUB  SO 
270  S0-S0+E"E 
275  El-E 

280  FOR  I-l  TO  N1 

285  B(I)-B(I)*1.001 

290  GOSUB  50 

295  B(I)-B(I)/1.001 

300  P(I)-(El-E)/( .001*8(1) ) 

305  NEXT  I 

310  FOR  I-l  TO  Nl 

315  G(I)-G(1)+E1*P(I) 

320  FOR  I 1-1  TO  Nl 

325  A(I,I1)-A(I,I1)  +  P(1)*P(I1) 

330  NEXT  II 
335  NEXT  I 
340  NEXT  J 
345  FOR  I-l  TO  Nl 
350  FOR  J-1  TO  Nl 

355  Q(I,J)-A(I,J)/(SQR( (A(I,I) )*(A(J, J) ) ) ) 

360  NEXT  J 

365  a(I)-0(I)/(8QR(A(l,I))) 

370  NEXT  1 

375  PRINT  I  PRINT  "S3E-" |S0)TAB( 22 )> "ITERATION  NO.-">T2; 
'i'AB(49)>*LAMBDA-"|L 

390  PRINT  "R-SQUARE  -" | 1- ( SO/ ( Z2- ( Zl *21/N) ) ) 

395  FOR  I-l  TO  Nl 
400  Q(I,I)«(J(I,I)*(1+].) 

405  NEXT  I 
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410  Q08UB  2250 

415  FOR  I-l  TO  N1 

420  P(I)-0 

425  FOR  J-1  TO  N1 

430  F(I)-P(I)+C(I,J)*a(J) 

435  NEXT  J 

440  P(T)-P(I)/(SOR(A(I,I))) 

445  NEXT  I 

450  IF  TKO  THEN  650 
455  IF  T2>-T1  THEN  500 
460  NN-0  t  I-l 
465  IF  I>N1  THEN  510 

470  IF  ABS<P(I>/B(I))<.0001  THEN  NN-NN+1 

475  IF  NN>-N1  THEN  490 

480  I-I4-1  I  OOTO  465 

490  PRINT  >  PRINT  "CONVERQBHCE" 

495  GOTO  620 

500  PRINT  "you  ARE  OUT  OF  ITERATIONS* 

505  GOTO  620 
510  FOR  I-l  TO  N1 
SIS  B<I)-B(I)'fP(I) 

520  NEXT  I 
525  81-0 

530  FOR  J-1  TO  N 
535  G08UB  50 
540  81-8148*8 
545  NEXT  J 

550  PRINT  "PARAMETERS" 

555  FOR  I-l  TO  N1 
560  PRINT  B(I)| 

570  NEXT  I 
575  PRINT 

580  IP  Sl>80  THEN  595 

585  L-L/10 

590  GOTO  230 

595  L-I.*10 

600  FOR  I-l  TO  Ml 

60S  B(I)-B(I)-P(n 

610  NEXT  I 

615  GOTO  395 

620  Tl— 1 

625  D-O 

630  GOTO  230 

650  PRINT 

655  V-80/(N-Nl) 

660  Vl-SQR(V) 

665  PRINT  *VAR-"fV|TAB<22))*SD-*;Vl|TAB(49>>*S8E-*)S0 

675  PRINT 

680  GOSUB  2275 

685  GOSUB  2300 

690  PRINT  "FINAL  PARAMETERS* ; TAB ( 22 "STD.  ERROR  OF  PARAM. 
TAB(49)|*CX)BPP.  OF  VAR." 

695  PRINT  " - •;TAB(22)»" - 

TAB(49)|" - " 

715  FOR  I-l  TO  Ml 
720  >SOR(A(I,I)) 

729  PRINT  TA8(3);B(I)>TAB(27>|DrrAB(S2)|D/B(I) 

730  NEXT  I 

735  PRINT  I  PRINT  "PRINT  VARIANCE-COVARIANCE  MATRIX 
(0-NO,l-yBS)*| 

740  INPUT  B7  i  PRINT 
745  IP  B7-0  THEM  765 
750  GOSUB  2325 
760  PRINT 

765  PRINT  "DO  VOU  NISM  PRINTOUT  OF  ESTIMATED  Y  FOR  EACH  X 
(0-N0,l-yBS)*| 

770  INPUT  H7  t  PRINT 
775  IF  H7-0  THEM  920 
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780  T-2 

785  PRINT  TAB(6)rx"»TAB(16)»“Y">TAB<23)»''ESTIMATED  Y"»TAB<36>> 
"a.E.  E8T.  Y">TAB(49)  j"RESIDUAI.'' 

790  PRINT  TAB(23)>"-(*;Tj"*  B.E. ) " jTAB( 36 ) | "+ ( " >T) S.E.)" 

792  PRINT  TAB(2)»'' - “iTABCiajj" - ";TAB(23); 

« - ''jTAB(36)  - '•»TAB(49)|" - " 

800  FOR  J«1  TO  N 
80S  GOSUB  50 
810  El-E 

815  FOR  I-l  TO  N1 

820  B(I)-B(I)*1.001 

825  GOSUB  50 

830  B(I)-B(I)/1.001 

835  P(I)-(E1-E)/(.001*B(I)) 

840  NEXT  I 
845  V1>0 

850  FOR  J2-1  TO  N1 
855  FOR  Jl-1  TO  N1 
860  Vl-Vl  *  A(.71,cr2)*P(Jl)*P(J2) 

865  NEXT  Jl 
870  NEXT  J2 
875  V2-SQR<V1) 

880  Ll-Y(a)-El-T*V2 
805  L2-L1+2*T*V2 

890  PRINT  TAB(2)iX(J,l))TAB(12)>Y(J>>TAB(23);Y(J)-El;TAO<36); 
V2>TAB(49) ;B1 

900  PRINT  TAB<23)|Ll|TAB(3e)|L2 
910  NEXT  J 
915  PRINT 

920  PRINT  "DO  YOU  WISH  TABLE  OF  RESIDUALS  ( O-NO r l-YES ) " » 

925  INPUT  H8 

930  IF  HS-O  THEN  970 

935  PRINT 

940  PRINT  "  RESIDUALS" 

945  PRINT 
950  FOR  J-1  TO  N 
955  GOSUB  50 
960  PRINT 
965  NEXT  J 

970  PRINT  «  PRINT  "PRINT  CORRELATION  MATRIX  (0-NO,1-YEB)" F 

975  INPUT  E8  i  IP  E0-O  THEN  1000 

980  PRINT 

985  GOSUB  2350 

1000  PRINT  I  INPUT  "PLOT  DATA  POINTS  (O-NO, 1“YES) " jDl  i  PRINT 
1005  IF  Dl-0  THEN  1500 

1500  INPUT  "PLOT  LEAST  SQUARES  PIT  ( O-NO, 1-YES ) " »D2  i  PRINT 
1505  IF  D2-0  THEN  100 
1995  END 

2000  REM  **  ZERO  MATRIX  G  •* 

2005  FOR  II-l  TO  N1 
2010  G(II)-0 
2015  NEXT  II 
2020  RETURN 

2025  REM  **  ZERO  MATRIX  A  ** 

2030  FOR  Il-l  TO  N1  «  FOR  JJ-l  TO  N1 
2035  A(II,J:)-0 
2040  NEXT  JJ  t  NEXT  II 
2045  RETURN 

2050  REM  **  ZERO  MATRIX  Q  ** 

2055  FOR  II-l  TO  N1  i  FOR  JJ-1  TO  N1 
2060  Q(II,JJ)-0 
2065  NEXT  JJ  l  NEXT  II 
2070  RETURN 

2075  REM  **  ZERO  MATRIX  C  ** 

2080  FOR  II-l  TO  N1  t  FOR  JJ-1  TO  N1 
2085  C(II,JJ>-0 
2090  NEXT  JJ  I  NEXT  II 
2095  RETURN 
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2100  REM  *  MATRIX  INVERSION  ROUTINE  * 

2105  FOR  JJ-1  TO  Nl 

2110  FOR  II-JJ  TO  Nl 

2115  IF  AA(II,JJ)<>0  THEN  2135 

2120  NEXT  II 

2125  PRINT  "SINGULAR  MATRIX" 

2130  STOP 

2135  FOR  KK-1  TO  Nl 
2140  S-AA(JJ,KK) 

2145  AA(JJ,KK)-AA(II,KK) 

2150  AA(II,KK)-S 
2155  S-C(JJ,KK) 

2160  C( JJ,KK)-C(1I,KK) 

2165  C(I1,KK)-S 
2170  NEXT  KK 
2175  TT-1/AA(JJ,JJ) 

2160  FOR  KK-1  TO  Nl 

2185  AA{JJ,KK)-TT*AA(JJ,KK) 

2190  C< JJ,KK)-TT*C<JJ,KK> 

2195  NEXT  KK 
2200  FOR  LL-1  TO  Nl 
2205  IF  LL-JJ  THEN  2235 
2210  TT— AA(LL,JJ) 

2215  FOR  KK-1  TO  Nl 

2220  AA(LL,KK)-AA(LL,KK)+TT*AA(JJ,KK) 

2225  C(LL,KK)-C(IiL,KK)+TT*C(JJ,KK) 

2230  NEXT  KK 
2235  NEXT  LL 
2240  NEXT  JJ 
2245  RETURN 

2250  REM  *  TRANSFER  MAT  Q  TO  MAT  AA  * 

2255  GOSUB  2075 

2260  FOR  II-l  TO  Nl  j  FOR  JJ-1  TO  Nl  »  AA( 11 , JJ >-Q ( 11 , JJ) 

2265  NEXT  JJ  t  C{II,I1>-1  i  NEXT  II 

2270  GOSUB  2100  t  RETURN 

2275  REM  *  TRANSFER  MAT  A  TO  MAT  AA  * 

2280  GOSUB  2075 

2285  FOR  ll-l  TO  Nl  «  FOR  JJ-1  TO  Nl  J  AA< II , JJ )»A( II , JJ) 
2290  NEXT  JJ  i  C(11,I1)-1  t  NEXT  11 
2295  GOSUB  2100  t  RETURN 

2300  REM  I  CONFIDENCE  MATRIX  TIMES  VARIANCE  * 

2305  FOR  ll-l  TO  Nl  t  FOR  JJ-1  TO  Nl 
2310  A(I1, JJ)-C(II,JJ)*V 
2315  NEXT  JJ,II 
2320  RETURN 

2325  REM  *  MATPRINT  A  * 

2330  FOR  II-l  TO  Nl  t  FOR  JJ-1  TO  Nl 
2335  E8-A(II,JJ)  :  PRINT  E8} 

2340  NEXT  JJ  :  PRINT  j  NEXT  11 
2345  RETURN 

2350  REM  *  MATPRINT  Q  * 

2355  FOR  II-l  TO  Nl  i  FOR  JJ-1  TO  Nl 
2360  E8-oai,JJ)  >  PRINT  B8; 

2365  NEXT  JJ  t  PRINT  i  NEXT  II 
2370  RETURN 

2400  REM  **  DATA  INPUT  ROUTINE  ** 

2405  INPUT  "ENTER  DATA  FILE  NAME">Z$ 

2410  OPEN  "I",1,Z$ 

2415  INPUT  #1,  N,N2 
2420  FOR  I-l  TO  N 
2425  INPUT  #1,  Y(U 
2430  FOR  J-1  TO  N2 
2435  INPUT  #1,  X(I,J) 

2440  NEXT  J,I 
2445  CLOSE  1 
2450  RETURN 
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APPENDIX  A 


PROGRAM  (3)  NONLIN/PRT 

5  LPRINT  "NONLINEAR  REGRESSION  ANALYSIS"  i  LPRINT 
31  LPRINT  "DATA  FILE  NAHE",,Z$ 

115  LPRINT  "MAX.  NO.  OF  ITERATIONS" ,T1 

125  LPRINT  "INITIAL  VALUE  OF  LAMBDA", L 

135  LPRINT  "NO.  OF  PARAMETERS* ,Nl 

141  LPRINT  I  LPRINT  "STARTING  PARAMETERS*  l  LPRINT 

150  PRINT  "B(*»I;")"j  t  INPUT  B(I)  «  LPRINT  "Bl " ) I j " ) " ,B( 1 )  i  NEXT  I 
175  LPRINT  :  LPRINT  i  LPRINT  "OBSERVATIONS  Y,X(l-5)"  l  LPRINT 
190  LPRINT  "<*»I;")"jTAB(7);y(I)»TAB(17);X(I,l)»TAB(27)jX(I,2)r 
TAB<37) jX(l,3) jTAB(47)jX(I,4>iTAB(57)fX(I,5) 

375  LPRINT  I  LPRINT  "SSE"* |S0)TAB( 22 )> "ITERATION  NO.-"|T2; 
TAB(49);"LAMBD’!-"»L 

390  LPRINT  -* >1- (S0/( Z2-( Z1*Z1/N ) ) ) 

490  LPRINT  I  I.rR/"T  "CONVERGENCE" 

500  LPRINT  "YO:'  >^‘.4E  OUT  OF  ITERATIONS* 

550  LPRINT  "PARAMETERS* 

560  LPRINT  B(I)r 

575  LPRINT 

650  LPRINT  I  LPRINT 

665  LPRINT  "VAR-* ;V)TAB(23) > *SD-"fVl}TAB( 49) ) "SSE-" >80 
690  LPRINT  I  LPRINT  l  LPRINT  *PI>>AL  PARAMETERS" >TAB( 22 ) > 

"STD.  ERROR  OF  PAHAM. * >TAB( 49) > "COEFF.  OP  VAR." 

695  LPRINT  " - "|TAB<22>>" - "» 

TAB(49)>" - - — * 

725  LPRINT  TAB< 3 ) >B< I ) >TAB(27 ) >D)TAB( 52 ) )D/B( I ) 

746  LPRINT  t  LPRINT  t  LPRINT  "VARIANCE-COVARIANCE  MATRIX"  i  LPRINT 
776  LPRINT  l  LPRINT  l  LPRINT  "PRINTOUT  OF  ESTIMATED  Y  FOR  EACH  X"  i 
LPRINT 

785  LPRINT  TAB( 6 )> "X" ;TAB(16 )) "Y")TAB( 23 )> "ESTIMATED  Y"| 
TAB(36)>"S.Ea  EST.  Y" >TAB( 49 )> "RESIDUAL" 

790  LPRINT  TAB<23))"-<*>T»""  S.E. )" »TAB( 36 >»"+(" (T) "*  S.E.)* 

792  LPRINT  TAB(2)j* - "jTAB(12)>" - ">TAB(23)> 

" - "»TAB(36>>" - ">TAB(49)»" - * 

890  LPRINT  TAB(2))X(J,l)>TAB(12)>y(J)>TAB(23}>y(J)-Sl>TAB(36)r 
V2>TAB(49)>E1 

900  LPRINT  TAB(23))L1)TAB(36)>L2 

940  LPRINT  I  LPRINT  i  LPRINT  "TABLE  OF  RESIDUALS*  t  LPRINT 
960  LPRINT  "("jJ>")",E 

981  LPRINT  I  LPRINT  i  LPRINT  "CORRELATION  MATRIX*  i  LPRINT 

1995  SYSTEM  "T*  t  END 

2335  xe-A(II,JJ)  t  LPRINT  EB) 

2340  NEXT  JJ  i  LPRINT  t  NEXT  II 
2360  ES-Q(II,JJ)  %  LPRINT  BS| 

2365  NEXT  JJ  i  LPRINT  t  NEXT  II 
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hPPGHDIX  A 

PROGRAM  (4)  AGAB/FUN  (1  PARAMETER) 

50  REM  **  FUNCTION  FOR  lal  ANTIGEN  ANTIBODV  BINDING  ** 

55  Pl-B(l)  i  TQ-X<J,1)*6.9E-6  i  TB-X( J,1 )*4.5E-6 
60  Al-Pl  «  Bl»l+Pl*Ta-Pl*TB  f  Cl— TB 
65  AB-(-Bl+8QRtBl*Bl-4*Al*Cl) )/(2*Al > 

70  F-P1*AB/(1+P1*AB) 

75  E-Y(J)-F 
80  RETURN 


PROGRAM  (4)  AGAB2/FON  (2  PARAMETERS) 

SO  reM  •*  FUNCTION  FOR  111  ANTIGEN  ANTIBODY  BINDING  ** 

55  Pl-B(l)  I  TG«X(J,l)*6.9B-6  t  TB«X(J,1)*B(2) 

60  Al-Pl  I  E1-1*'P1*TG-P1*TB  I  Cl— TB 
65  AB-(-alf8QR(Bl*Bl-4*Al*Cl) )/(2*Al» 

70  F-P1*AB/<1+P1*AB) 

75  E-Y(J)-F 
SO  RETURN 


PROGRAM  <4)  AB2AG/FUN  (3  PARAMETERS) 

50  REM  *  FUNCTION  FOR  2  DISTINCT  SPECIES  OF  ANTIBODY,  DIFFERENT  K'S 
52  Pl-Bd)  1  P2-B«2)  I  TB-X(J,l)*4.5E-6  «  T(l-X  ( J,  1 )  *0 . 9E-6 
54  BL>0  t  LQ-0  t  TL-TG 
56  R«SIN(B(3))*SIN(B(3)) 

58  0«..5’*(BL+TI,) 

60  B1-R*TB/(1+P1*G) 

62  B2-(1-R)*TB/(1+P2*G) 

64  TC-G*(1+P1*B1+P2*B2) 

66  IF  TO TO  THEN  TD-G  ELSE  BL-O 
68  IP  ABS( (LG-G)/a)<lE-6  THEN  72 
70  LQ-G  I  GOTO  58 

72  F-(TG-0)/TQ  J  E-Y(J)-F  i  RETURN 


APPENDIX  B 


EXAMPLE  <1)  -  FIXED  TOTAL  ANTIGEN  AND  ANTIBODY  CONCENTRATION 


50  REi*  **  FUNCTION  FOR  Itl  ANTIGEN  ANTIBODY  BINDING 
55  Pl>B(l)  «  TO«X<J,L)*6.9E-6  i  TB-X( J,1 )*4.5B-6 
60  Al-Pl  «  Bl-l+Pl*Ta-Pl*TB  t  Cl"-TB 
65  AB- ( -Bl+SQR( BL*B1-4*A1*C1 ) )/ < 2*A1 ) 

70  P-Pl*AB/(l+Pl*AB) 


75  F.-Y(J)-F 
60  RETURN 


<1 


NONLINEAR  REGRESSION  ANALYSIS 


DATA  FILE  NAME 
MAX.  NO.  OF  ITERATIONS 
INITIAL  VALUE  OF  LAMBDA 
NO.  OF  PARAMETERS 


AGAB2 01/DAT 
15 
1 
1 


STARTING  PARAMETERS 


B(  1  )  lE+06 


OBSERVATIONS  Y, 

,X(l-5) 

(  1  ) 

.336 

1 

0 

(  2  ) 

.271 

.5 

0 

(  3  ) 

.247 

.25 

0 

(  4  ) 

.232 

.2 

0 

(  5  ) 

.231 

.166667 

0 

(  6  ) 

.21 

.125 

0 

(  7  ) 

.196 

.1 

0 

(  8  ) 

.147 

.05 

0 

(  9  ) 

.122 

.0333333 

0 

(  10  ) 

.108 

.025 

0 

(  11  ) 

.093 

.02 

0 

SSE-  .0788429  ITERATION  NO.-  1 

R-SQUARE  —.391782 

PARAMETERS 

723492 


SSE-  ,0440415  ITERATION  NO.-  2 

r-SQUARE  -  .22255 

PARAMETERS 

534537 


SSE-  .032425  ITERATION  NO.-  3 

R>SQUARB  -  .427612 

PARAMETERS 

517085 


SSE-  .0322849  ITERATION  NO.-  4 

R-SQUARB  -  .430085 

PARAMETERS 

514760 


BSE-  .0322822  ITERATION  NO.-  5 

R-SQUARB  -  .430132 

PARAMETERS 

914392 

SSE-  .0322822  ITERATION  NO.-  6 

R-SQUARB  -  .430133 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


LAMBDA-  1 


LAMBDA-  .1 


LAMBDA-  .01 


LAMBDA-  lE-03 


LAMBDA-  lE-04 


LAMBDA-  lE-05 


PARAMETER) 


CONVERGENCE 

SSE-  .0322822 
R-SQUARB  -  .430133 


ITERATION  NO.-  7 


LAMBDA-  0 


APPENDIX  B 


EXAMPLE  (1)  -  FIXED  TOTAL  ANTIGEN  AND  ANTIBODY  CONCENTRATION  -  ( 


VAR-  3.22822E-03 

SD-  .0568174 

SSB-  .0322822 

FINAL  PARAMETERS 

STD.  ERROR  OF  PARAM. 

COEFF.  OF  VAR, 

514392 

87748.9 

.170588 

VARIANCE-COVARIANCE  MATRIX 
7.69987B+09 


PRINTOUT  OP  ESTIMATED  Y  FOR  EACH  X 


X 

Y 

ESTIMATED  Y 
-<  2  *  S.E. 

S.E.  BST.  Y 
>  +(  2  *  S.E. 

RESIDUAL 

) 

1 

.336 

.435135 

.0196493 

-.0991355 

.395837 

.474434 

.5 

.271 

.349467 

.0221201 

-.0784667 

.305227 

.393707 

.25 

.247 

.258758 

;0219879 

-.0117584 

.214783 

.302734 

.2 

.232 

.23041 

(.0212914 

1.58979E-03 

.187827 

.272993 

.166667 

.231 

,208058 

.0204907 

.0229425 

.167076  . 

.249039 

.125 

.21 

.174782 

.0186867 

.0352179 

,137009 

.212556 

.1 

.196 

,151015 

.0174022 

.0449853 

.11621 

.185819 

.05 

.147 

.0906243 

.012251 

.0563757 

.0661224 

.115126 

.0333333 

.122 

.0649603 

9.38744B-03 

.0570397 

.0461854 

.083/352 

.025 

.108 

.0506688 

7.59664E-03 

.0573312 

.0354753 

.0656621 

.02 

.093 

.0415453 

6.37967E-03 

.0514547 

.0287859 

.0543046 

TABLE  OF 

RESIDUALS 

(  1  ) 

-  ^ 

0991355 

(  2  ) 

- 

0704667 

(  3  ) 

- 

0117584 

<  4  ) 

,.50979E-O3 

(  5  ) 

0229425 

(  6  ) 

0352179 

(  7  ) 

0449853 

(  8  ) 

0563757 

(  9  ) 

0570397 

(  10  ) 

0573312 

(  11  ) 

0514547 

I  PARAMETER) 


CORRELATION  MATRIX 
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EXAMPLE  (2)  -  VARIABLE  TOTAL  ANTIBODY  CONCENTRATION  -  (2  PARAMETERS) 

50  REM  **  FUNCTION  FOR  111  ANTIGEN  ANTIBODY  BINDING  ** 

55  P1-B<1)  t  TG«X(J,I)*6.9E-6  :  TB««X{J,1)*B(2) 

60  Al-Pl  I  Bl-1+P1*TG-P1*TB  i  Cl—TB 
65  AB-(-Bl+SQR(Bl*Bl-4*Al*Cl) )/(2*Al) 

70  F-n*AB/(l+Pl*AB) 

75  E-Y(J)-F 
80  RETURN 


NONLINEAR  REGRESSION  ANALYSIS 


DATA  FILE  NAME 
MAX.  NO.  OF  ITERATIONS 
INITIAL  VALUE  OF  LAMBDA 
NO.  OF  PARAMETERS 


AGAB20 1/DAT 
15 
1 
2 


STARTING  PARAMETERS 


B(  1  )  lE-fOe 

B(  2  )  4.5E-06 


OBSERVATIONS  Y,X(.l-5) 


(  1  ) 

.336 

1 

0 

(  2  ) 

.271 

.5 

0 

(  3  ) 

.247 

.25 

0 

(  4  ) 

.232 

.2 

0 

(  5  ) 

.231 

.166667 

0 

(  6  ) 

.21 

.125 

0 

(  7  ) 

.196 

.1 

0 

(  e  ) 

.147 

.05 

0 

(  9  ) 

.122 

.0333333 

0 

(  10  ) 

.108 

.025 

0 

(  11  ) 

.093 

.02 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


SSE>  .0788425  ITERATION  HO.-  1  LAMBDA-  1. 

R-8QUARE  —.391782 

PARAMETERS 

856896  3.8962E-06 


SSE-  .0271986  ITERATION  HO.-  2  LAMBDA-  .1 

R-SQUARE  -  .519872 

PARAMETERS 

1.13332E+0e  2.B4243E-06 


SSE-  .0103209  ITERATION  NO.-  3  LAMBDA-  .01 

R-SQUARE  -  .817809 

PARAMETERS 

1.93685B4-06  2.22722E-06 


SSE-  7.65994B-03  ITERATION  NO.-  4  LAMBDA-  lE-03 

R-SQUARB  -  .864782 

PARAMETERS 

2.S3919B406  2.292e8E-06 


SSE-  1.6S261E-03  ITERATION  NO.-  5  LAMBDA-  lE-04 

R-SQUARE  -  .970827 

PARAMETERS 

2.6S473B-t-Q6  2.2902SE-06 

SSE-  1.57879E-03  ITERATION  HO.-  6  LAMBDA-  lE-05 

R-SQUARB  -  .97213 

PARAMETERS 

2.65  934E-t-06  2.28993E-06 


APPENDIX  B 


EXAMPLE  (2)  -  VARIABLE  TOTAL  ANTIBODY  CONCENTRATION  -  (2  PARAMETERS) 

SSE-  1.57867E-03  ITERATION  NO.-  7  LAMBDA-  lE-06 

R- SQUARE  -  .972132 

PARAMETERS 

2.e6445E4-06  2.2ea95B-06 

SSE-  1.57865E-03  ITERATION  NO.-  8  LAMBDA-  lE-07 

R-SQUARE  -  .972133 

PARAMETERS 

2.67128E4-06  2.2874E-06 

SSE-  1.5786E-03  ITERATION  NO.-  9  LAMBDA-  lE-08 

R-SQUARE  -  .972134 

PARAMETERS 

2.66406E4-06  2.28e76B-06 

SSE-  1.578S7E-03  ITERATION  NO.-  10  LAMBDA-  lB-09 

R-SQUARE  -  .972134 

PARAMETERS 

2.65B47E4-06  2.2900SE-Q6 

PARAMETERS 

2,6S847E-f06  2.29005E-06 

PARAMETERS 

2.65847E406  2.29005E-0e 

PARAMETERS 

2,6S847E-f06  2.2900SE-a6 

PARAMETERS 

2.65847E-t'Q6  2.2900SE-06 

PARAMETERS 

2.65a47E+06  2.29005E-06 

PARAMETERS 

2.6585E't-06  2.29004E-06 

PARAMETERS 

2.658aiE'fC6  2.28996E-06 

PARAMETERS 

2,66064E-t-06  2.28949E-06 

PARAMETERS 

2.66324E-t-06  2.2e887E-06 

CONVERGENCE 

SSE-  1.37857E-03  ITERATION  NO.-  11  LAMBDA-  0 

R-SQUARE  -  .972134 


VAR-  1.75397E-04  SD-  .0132437  SSE-  1.57857B-03 


FINAL  PARAMETERS  STD.  ERROR  OF  PARAM.  COBPF.  OP  VAR. 


2.66406E-t-06  322641  .121109 

2.28e76B-06  8.48315E-08  .0370644 


VARIANCE-COVARIANCE  MATRIX 

1.04097E'fll  -.0229482 
-.0229482  7.19637B-15 


PRINTOUT  OP  ESTIMATED  Y  FOR  EACH  X 

X  Y  ESTIMATED  Y  8.E.  BST.  Y  RESIDUAL 

-(  2  *  3.E. )  +(  2  *  S.E. ) 


1  .336  .307643  9.02226E-03  .0284566 

.289499  .325588 
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EXAMPLE  (2)  -  VARIABLE  TOTAL  ANTIBODY  CONCENTRATION  -  (2  PARAMETERS) 


.,5 

.271 

.287748 

.274019 

6.86438B-03 

.301477 

-.0167481 

.25 

.247 

.256595 

.246823 

4.88633E-03 

.266368 

-9.59546E~03 

.2 

.232 

.243942 

.234855 

4.5436E-03 

.253029 

-.011942 

.166667 

.231 

.23271 

.223824 

4.44295E-03 

.241596 

-1.70967E-03 

.125 

.21 

.213538 

.204382 

4.57605E-03 

.222694 

-3.53788E-03 

.1 

.196 

.197674 

.187976 

4.848648-03 

.207371 

-1.67352E-03 

.05 

.147 

.145878 

.13470] 

5.58837E-03 

.157055 

1.1219iE-03 

.0333333 

.122 

.1165 

.105429 

5.53572E-03 

.127572 

5.49953E-02 

.025 

.108 

.0972603 

.0867616 

5.2493SE-03 

.107759 

.0107398 

.02 

.093 

.0835914 

.0737701 

4.91065E-03 

.0934127 

9.4O858E-03 

TABLE  OF 

RESIDUALS 

(  1  ) 

.0284566 

(  2  ) 

-.0167481 

(  3  ) 

-9.59546E-03 

(  4  ) 

-.011942 

<  5  ) 

-1.70967E-03 

(  6  ) 

-3.53788E-03 

<  V  ) 

-1.67352E-03 

(  a  ) 

1.12191E-03 

(  9  ) 

5.499S3E-03 

<  10  ) 

.0107398 

(  11  ) 

9.40858E-03 

CORRELATION  MATRIX 

.999997  .638443 

.838443  1 
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EXAMPLE  (3)  -  2  SPECIES  ANTIBODY,  2  BXNDINQ  CONSTANTS  -  (3  PARAMETERS) 

50  REM  **  FUNCTION  FOR  2  DISTINCT  SPECIES  OF  ANTIBODY,  DIFFERENT  K'S  ** 
52  Pl-B(l)  I  P2-B(2>  I  TG-X(J,l)*6.9E-6  i  TB-X( J,l) *4.5E-6 
54  BL-0  I  LG«0  I  TL-TQ 
56  R-SIN(B(3>>*SIN(B(3> ) 

58  G-.5*<BL+TL1 
60  Bl-R*TB/(1+Pl*a) 

62  B2-(1-R)*TB/(1+P2*G) 

64  TC-Q*(1+P1*B1+P2*B2) 

66  IP  TC>TO  THEN  TL»G  ELSE  BL»‘G 
68  IF  ABS( (LG-a)/G><lB-6  THEN  72 
70  LG-G  I  GOTO  58 

72  F-(TG-G)/TG  t  E-Y(J)-P  t  RETURN 


NONLINEAR  REGRESSION  ANALYSIS 

DATA  PILE  NAME  AGAB20I/DAT 

MAX.  NO.  OF  ITERATIONS  15 

INITIAL  VALUE  OF  LAMBDA  1 

NO.  OF  PARAMETERS  3 

STARTING  PARAMETERS 

B(  1  )  1E406 

Bt  2  )  10000 

B(  3  )  .4 


OBSERVATIONS  Y,X(l-5) 


(  1  ) 

.336 

1 

0 

0 

0 

0 

(  2  ) 

.271 

.5 

0 

0 

0 

0 

(  3  ) 

.247 

.25 

0 

0 

0 

0 

(  4  ) 

.232 

.2 

0 

0 

0 

0 

(  5  ) 

.231 

.166667 

0 

0 

0 

0 

(  6  ) 

.21 

.125 

0 

0 

0 

0 

(  7  ) 

.196 

.1 

0 

0 

0 

0 

(  8  ) 

.147 

.05 

0 

0 

0 

0 

(  9  ) 

.122 

.0333333 

0 

0 

0 

0 

(  10  ) 

.108 

.025 

0 

0 

0 

0 

(  11  ) 

.093 

.02 

0 

0 

0 

0 

SSE-  . 

259762 

ITERATION 

NO.-  1 

LAMBDA- 

1 

R-SQUARK  — 3.5R55 
PARAMETERS 

3.42768E-t-06  37845.6  .589247 

SSE-  .0276888  ITERATION  NO.-  2  LAMBDA-  .1 

R-SQUARE  -  .511219 

PARAMETERS 

4.98381E4-06  62957.1  .649253 

SSE-  6.29858B-04  ITERATION  NO.-  3  LAMBDA-  .01 

R-SOUARE  -  .988882 

PARAMETERS 

4.47224E4-06  62553.9  .670344 

SSE-  3.31247E-04  ITERATION  NO.-  4  LAMBDA-  lE-03 

R-SQUARE  -  .994153 

PARAMETERS 

4.42239B-I-06  60377.4  .674733 

SSE-  3.26629E-04  ITERATION  NO.-  5  LAMBDA-  lE-04 

R-SQUARE  -  .994234 

PARAMETERS 


EXAMPLE!  O)  -  2  SPECIES  ANTIBODY,  2  BINDING  CONSTANTS 


APPENDIX  B 
(3  PARAMETERS) 


4.43023E+Oe  60681.1 

PARAMETERS 

.6743 

4.4299E+06  60673.7 

PARAMETERS 

.674314 

4,42''52E-t'06  60619.5 

PARAMETERS 

.674418 

4.42257E+06  60491.7 

.674649 

SSE-  3.26614B'04 

R  'SQUARE  *  .9  9423  4 

ITERATION  NO 

6 

LAMBDA-  .01 

PARAMETERS 

4.40268E')-06  59897.8 

PARAMETERS 

.675587 

4.4181?E-f06  60266.9 

PARAMETERS 

.674898 

4  12229E4-06  60437.9 

.674663 

SSE-  3.26613E-04 
R-SQUr.RB  -  .994234 
PARAMETERS 

ITERATION  NO 

.-  7 

LAMBDA-  .1 

4.4?lS6E-t-06  60452.2 

PARAMETERS 

.67469 

4.42223E-t-06  60447.9 

.674671 

'T'.'NVERGENCE 

SSE-  3.266iaE-04 
R-SQUARE  -  .994234 

ITERATION  NO 

B 

LAMBDA-  0 

VAR-  4.08266E‘-05 

SD-  6.3895BE 

-03 

SSE-  3.26613E‘ 

FINAL  PARAMETERS 

STD.  ERROR  OP  PARAM. 

COEFF.  OF  VAR 

4.42229E'l-0e 

550569 

.124499 

60437.9 

12239.4 

.202512 

.674663 

.0223403 

.0331133 

v^ARIANCE-COVARIANCE  MATRIX 

3.03127E-(11  5.1076E+09  -11292.8 

S.J076E+09  1.49804E+08  -248.384 

-11292.8  -248.384  4. 

9909E-04 

PRINTOUT  OF  ESTIMATED 

Y  FOR  EACH  X 

X  Y 

ESTIMATED  Y 

S.E.  EST.  Y 

RESIDUAL 

-(  2  *  S.E.  ) 

+(  2  *  S.E. 

) 

1  .336 

.329505 

5.84322B-03 

e.49476E-03 

.317819 

.341192 

.5  .271 

.284716 

3.2307B-03 

-.013718 

.278255 

.291177 

.25  .247 

.245774 

3.024eE-03 

1.22634E-03 

.239724 

.251823 

.2  .232 

.233536 

2.93133E-03 

-1.53597E-03 

.227673 

.239399 

.166667  .231 

.223414 

2.7993E-03 

7.58627E--03 

.217815 

.229012 

.125  .21 

.207019 

2.59741B-03 

2.9eOBlE-03 

.201824 

.212214 

.1  .196 

.193848 

2.49761E-03 

2.15201E-03 

.188853 

.198843 

E 

.-.■i 

1-7 
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EXAMPLE  (3)  -  2  SPECIES  ANTIBODY,  2  BINDING  CONSTANTS  -  (3  PARAMETERS) 


.05  .147 
.0333333  .122 
.025  .108 
.02  .093 


.150552 

.145066 

.12463 

.118288 

.106789 

.100232 

.0936081 

.0069725 


2.74308E-03  -3.55201E-03 
.156038 

3.17106E-03  -2.6303E-03 
.130972 

a.27849E-03  1.21116E-03 

,113346 

3.31783B-03  -6.08131E-04 
.100244 


TABLE  OF  RESIDUALS 

6.49476E-03 
.013716 
1.22634E-03 
1.53597E-03 
7.58627E-03 
2.9eoeiE-03 

2.15201E-03 
3.55201E-03 
2.6303E-03 
1.21116E-03 
'6.08131E-04 


CORRELATION  MATRIX 

.999996  .459018  .841831 

.459018  1  .822007 

,841831  .822007  1 


(  1  ) 

(  2  ) 

(  3  ) 

(  4  ) 

(  5  ) 

(  6  ) 

(  7  ) 

(  8  ) 

(  9  ) 

(  10  ) 
{  11  ) 


R-  .390168 


(1-R)«  .609832 


APPENDIX  C 


DERIVATION  OF  111  ANTIOEN-ANTIBODY  BIKDINO  FUNCTION 
Definitions i 

G>free  antigen  cono.  B>free  antibody  cone. 

TG>total  antigen  cone.  TB^total  antibody  oonc. 

BQ»baund  antlgan-antlbody  oomplex  K»blnding  constant 

Equilibrium  expressiont 

K-Ba/(B*G)  or  BG-K*B*G  (1) 

Masn  balance  equationai 

TB-B-t'BG  (2)  and  TG-G+BQ  (3) 


Using  equations  (1)  and  (3)  we  geti 

T0-O+K*B*a  or  TG-a*tl+K*B)  or  G-T0/(1+K*B)  (4) 
Using  equations  (2)  and  (3)  we  geti 

TB-B+K*B*G  and  using  (4)i  TB-B+(K*B*TG/(1+K*E) )  (5) 

Using  equations  (2)  and  (5)  we  geti 

BG-TB-B  or  Ba-B+<K*B*TG/(1+K*B> )-B  or  BG"K*B*Ta/( 1+K*B) 
Dividing  both  sides  by  TO  we  geti 
(BG/TG)-K*B/(1+K*B)  (6) 


(BG/TG)  Is  the  fraction  antigen  bound  or  our  measured  Y 
value.  If  we  let  K*B  equal  X,  we  see  that  the  lil  binding 
equation  reduces  to  the  hyperbolic  form  Y^X/d-fX),  where 
Yafraction  antigen  bound  and  Xathe  binding  constant  times  the 
free  antibody  concentration. 

In  order  to  solve  equation  (6)  for  fraction  antigen  bounds 
we  must  first  know  the  free  antibody  concentration  (B).  We  can 
use  equation  (5)  which  contains  only  known  quantities  and  B«  and 
reduce  it  to  a  second  order  polynomial  In  B. 

Multiplying  each  term  of  equation  (3)  by  (l'fK*B)  we  geti 

TB*(1+K*B)-B*{1+K*B>+K*B*T0  or  TB+K*B*TB-B+K*B*B+K*D*TO 

Rearranging  and  collecting  terms  we  geti 

K*(B*B)+B+K*B*TO-K*B*TB-TB-0  or  K*(B*B)+(l+K*(TO-TB) >*B-TB-0 

This  la  a  quadratic  equation  in  B  where  the  terns  atei 
a-M  b«l-»-R*(TO-TB) )  and  c— TB 

It  Is  now  possible  to  use  the  quadratic  aquation 
((-bH-BQR(b*b->4*a*o))/2*a)  to  solve  for  B  unlng  these  terms. 

These  are  the  same  terms  used  in  function  AGAB/FUN  In  Appendix  A 
and  In  examples  1  and  2  Appendix  B.  Once  B  Is  known  then 
aquation  (C)  can  be  used  to  oaloulate  the  fraction  antigen  bound 
and  the  error  tarm  to  complete  the  function  subroutine.  This 
sane  proceduru  can  be  used  to  solve  other  equilibrium  models 

Provided  a  means  for  finding  the  root  of  the  derived  polynomial 
a  available. 


APPENDIX  C 


MODEL  FOR  ANTIQEN-ANTIBODY  BINDING  ASSUMING  TWO 
DISTINCT  SPECIES  OF  ANTIBODY  WITH 
DIFFERENT  BINDING  CONSTANTS  FOR  ANTIGEN 


Daflnitlonat 

G-frae  anblgen  cone. 
TBl-total  cone,  antibody  1 
TB2ototal  cone,  antibody  2 
TB~total  antibody  cone. 
TG-total  antlgan  cone, 
GBlocomplex  with  Bl 


Bl"fraa  antibody  1  cone. 
B2'«£raa  antibody  2  cone. 
R-fractlon  of  TB  aa  TBl 
Kl-blndlng  constant  1 
K2°blndlng  constant  2 
GBiMcomplax  with  B2 


Equilibrium  axpraaalonat 
aBl-Kl*a*Bl  and  aB2-K2*0*B2 
Mass  balanca  aquations) 


TBl-Bl+GBl 

and 

TB2-B2'faB2 

and 

TQ-Q+GB1+GB2 

Partition 

functional 

TB-TB1+Ta2 

and 

TB1-R*TB 

and 

TB2-(1-R)*TB 

Using  tha  equilibrium  and  mass  balanca  equations  we  geti 


Bl-TB1/(1+Kl*a) 

or 

B1-R*TB/(1+K1*G)  (1) 

B2-TH2/'(1+K2*a) 

or 

B2-(l-R)*TB/a+K2*G) 

(2) 

Ta-G+Kl*a*Bl+K2*G*B2 

or 

TG-G*(1+K1*B1+K2*B2) 

(3) 

If  wa  ware  bo  solva  thasa  aquations  in  barms  of  Q  wa  would 
obtain  a  third  order  polynomial  or  cubic  aquation  in  G.  since 
there  is  no  simple  technique  for  obtaining  the  root  of  this 
equationi  there  is  no  point  in  solving  in  terms  of  G.  Instead  we 
use  a  numariujl  technique  for  solution  which  depends  on  the 
equations  being  positive  definite  functions. 

Equations  (l)-(3)  above  are  structured  so  that  given  G  we 
can  solve  for  TG.  TQ  is  known  however  and  is  the  key  to  knowing 
when  we  have  selected  the  correct  G.  First  wa  sat  a  lower  and  an 
upper  limit  to  G.  We  set  the  lower  limit  to  0  and  the  upper 
limit  to  TG.  We  then  guess  that  the  value  for  G  is  half  way 
between  0  and  TQ.  This  first  value  of  Q  is  used  to  calculate  a 
value  for  TG  and  compare  it  to  the  real  value.  If  the  calculated 
value  is  too  high  then  B  was  too  high  and  we  make  tha  first  guess 
on  G  the  new  upper  limit.  If  it  is  too  low  wo  make  the  first 
guess  on  G  the  new  lower  limit.  We  make  a  new  guess  half  way 
between  the  new  upper  and  lower  limits  and  try  again  until  the 
value  for  G  converges  to  the  correct  value.  This  usually  happens 
in  7-10  passes.  Once  we  know  the  value  for  0  we  can  calculate 
the  fraction  bound  ( F>(TG-a)/Ta) . 
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Th«  following  BASIC  oodo  perfocma  the  numarioal  aolution  of 
a  function  ualng  thla  binary  aaaroh  mathodi 

50  REM  **  PUNCTIOH  FOR  2  DIFFERENT  SPECIES  OF  ANTIBODY  WITH 
DIFFERENT  K'B  -  Pl-Kl  AND  P2-K2  ** 

52  Pl-B(l)  <  P2>B(2>  I  TB-X(J,l)*4.5B-6  t  TQ-XC J,I)*6.9E-6 
54  BI.X0  I  IiO«0  I  TL-TO 
56  R«SIN(B(3))*8IN(B(3)) 

50  0-.5*<BL+TL) 

60  Bl-R*TB/(1+Pl*a) 

62  B2-<1-H)*TB/(1+P2*a) 

64  TC-a*(l+Pl*Bl+P2*B2) 

66  IF  XC>TQ  THEN  TL-Q  ELSE  BL-G 
68  IF  AB8( (LQ-a)/a)<lB-6  THEN  72 
70  LQ-Q  I  GOTO  SB 

72  F«(Ta-0)/TG  I  E>y(J)-F  i  RETURN 

Tha  varlablaa  art  daflnad  aa  followai 

Pl-Kl r  P2-K2)  R-  fraction  total  antibody  aa  apaclaa  Bl. 

B(l)i  B(2)  and  B(3>  ara  tha  paramatar  varlablaa  B'a 

manlpulatad  by  tha  main  program  rapraaanting  XI,  K2, 
and  R  raapaotlvaly. 

TB  and  TO  ara  tha  total  antibody  and  total  antlgan 
ooncantratlona. 

X(J,1>  la  tha  X  data  varlabla,  J  la  tha  numbar  of  tha  data 
point  and  1  maana  tha  flrat  Indapandant  varlabla  of  a 
poaalbla  5. 

Q  la  tha  fraa  antlgan  conoantrablon. 

Bt  la  tha  bottom  limit  for  G,  TL  la  tha  top  limit  for  G  and 
LG  la  tha  laat  valua  for  G. 

Bl  and  B2  ara  tha  fraa  ooncantratlona  of  antibodlaa  1  and 
antibody  2  raapactlvaly. 

TC  la  tha  calculatad  concantration  for  TG,  Ta(calo)4 
Y(J)  la  tha  Y  valua  for  tha  Jth  data  point. 

Lina  56  provldaa  a  naw  guaaa  for  G  half  way  batwaan  tha 
uppar  and  tha  lowar  limits. 

Llnaa  60  to  64  calculata  a  valua  for  TOloalc)  ualng  tha  G 
from  llna  58. 

Lina  66  comparaa  tha  TGloalo)  with  TO  and  sata  tha  naw  uppar 
or  lowar  limit  dapandlng  on  tha  outooma. 

Lina  66  taata  for  convarganca  and  althar  branohaa  out  of  tha 
loop  or  cauaaa  anothar  paaa. 

Lina  70  rasats  tha  last  Q  valua  to  tha  currant  valua. 
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VAHIABLE  CROSS  REFERENCE  INDEX  FOR  "MONLIN/PRT* 

A  25(  *325(2  355(3  36S(  440(  720(  860(  *2035(  22e5(  *2310(  233S( 

AA  25(  2113(  2140(  *2145(2  *2150(  2175(  *2185(2  2210(  *2220(3  *226'0(  *22eS( 

B  20(  *150(2  *285(2  *295(2  300(  470(  *515(2  S60(  *605(2  725(2  *820(2  *830(2 
835  ( 

C  25(  430(  *2085(  2155(  *2160(2  *2165(  *2190(2  *2225(3  *2265(  *2290(  2310( 

D  *720  725/2 
D1  *1000  1005 
D2  *1500  1505 

E  270/2  275  300  540/2  810  835  960 
El  *275  300  315  *810  835  880  890/2 
B5  *160  170 
B6  *100  105 
E7  *740  745 

E8  *975/2  *2335/2  *2360/2 
Q  20(  *315(2  *365(2  43Q(  *20I0( 

H7  *770  775 
H8  *925  930 

I  15  *145  150/5  *180  190/7  195  *205  210  215/2  220  *280  285/2  295/2  300/2  305 
*310  315/3  325/3  335  *345  355/4  365/4  370  *395  400/4  405  *415  420  430/3 
440/4  445  *460  465  470/2  *480/2  *510  515/3  520  *555  560  570  *600  605/3  610 
*715  720/2  725/2  730  *815  820/2  830/2  835/2  840  *2420  2425  2435  2440 

II  320  325/3  330 

II  *2005  2010  2015  *2030  2035  2040  *2055  2060  2065  *2080  2085  2090  *2110  2115 
2120  2145  2150  2160  2165  *2260/3  2265/3  *2285/3  2290/3  *2305  2310/2  2315 
*2330  2335  2340  *2355  2360  2365 

J  15  *260  340  *350  355/4  360  *425  430/2  435  *530  545  *800  880  890/3  910  *950 
960  965  *2430  2435  2440 

31  *855  860/2  865 

32  *850  860/2  870 

33  *2030  2035  2040  *2055  2060  2065  *2080  2085  2090  *2105  2110  2115  2140  2145 
2155  2160  2175/2  2185/2  2190/2  2205  2210  2220  2225  2240  *2260/3  2265  *2285/3 
2290  *2305  2310/2  2315  *2330  2335  2340  *2355  2360  2365 

K  15 

KK  *2135  2140  2145/2  2150  2155  2160/2  2165  2170  *2180  2185/2  2190/2  2195  *2215 
2220/3  2225/3  2230 

I,  *120  125  375  400  *585/2  *595/2  *625 
Ll  *880  885  900 
L2  *885  900 

LL  *2200  2205  2210  2220/2  2225/2  2235 
N  180  205  260  390  530  655  800  950  *2415  2420 

Ml  *130  135  145  280  310  320  345  350  395  415  425  465  475  510  555  600  655  715  815 
650  655  2005  2030/2  2055/2  2080/2  2105  2110  2135  218C  2200  2215  2260/2 
2285/2  2305/2  2330/2  2355/2 
N2  *2415  2430 
NN  *460  *47C/2  475 

P  20(  *300(  315(  325(2  *420(  *430(2  *440(2  470 (  515(  605 (  *a35(  860(2 
Q  25(  *355 (  *400(2  *20e0(  2260(  2360( 

S  *2140  2150  *2155  2165 

80  *235  *270/2  375  390  580  655  665 

81  *’>25  *540/2  580 

T  */80  790/2  880  885 
T1  *110  115  450  455  *620 
T2  *225  *23(J/2  375  455 
TT  *2175  2185  2190  *2210  2220  2225 

V  *655  660  $65  2310 

VI  *660  665  *845  *860/2  875 

V2  *875  880  865  890 

X  20(  190(5  890(  *2435( 

X  20(  190(  210(  215(2  B80(  890(2  *2425( 

Z  31/8  *2405/8  2410/8 
Z1  *200  *210/2  390/2 
Z2  *200  *215/2  390 
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